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ABSTRACT 

Based on observations of the most abundant sulfur-bearing molecules (H2S, 
CS, NS, SO, H 2 CS, OCS, and S0 2 ) carried out with the IRAM 30m telescope 1 
and SEST 2 , we present the first analysis of the sulfur chemistry of an extragalac- 
tic source, the nuclear region of the starburst galaxy NGC 253. This is the first 
time that H 2 S and, tentatively, H 2 CS are detected towards the nucleus of a star- 
burst galaxy Source averaged fractional abundances of these molecules are a few 
10~ 9 , except for CS and OCS which are more abundant (10~ 8 ). Sulfur isotopic 
ratios, 32 S/ 34 S~8±2 and 34 S/ 33 S>9, are measured through observations of 13 CS, 



1 Based on observations carried out with the IRAM 30m telescope. IRAM is supported by INSU/CNRS 
(France), MPG (Germany) and IGN (Spain). 

2 Based on observations with the Swedish/ESO submillimetre telescope (SEST) at the European Southern 
Observatory, Chile. ESQ N°70.B-0022. 
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C 34 S, and C 33 S. A comparison with the observed relative abundances towards 
different prototypical Galactic sources suggests that the chemical composition of 
NGC253 is similar to that found towards the molecular clouds complexes like 
SgrB2 in the nuclear region of the Milky Way. The large overabundance of OCS 
compared to the predictions of time-dependent sulfur chemistry models supports 
the idea that OCS is likely injected into the gas phase from the grain mantles by 
low velocity shocks. 

Subject headings: galaxies: abundances — galaxies: active — galaxies: individ- 
ual (NGC 253) — galaxies: ISM — galaxies: nuclei — galaxies: starburst 

1. Introduction 

The formation and evolution of stars strongly affects the chemistry of their molecular 
environment. In particular, sulfur-bearing molecules show enhanced abundances in the hot 
molecular cores, the condensations associated with massive protostars. The most accepted 
idea is that during the collapse of a star forming cloud, a large fraction of the available 
sulfur freezes out by accretion onto grain mantles and hydrogenates in the form of H2S. The 
high temperatures and turbulent motions in these objects produced by the newly formed 
stars induce the evaporation and disruption of grain mantles, releasing the H 2 S into the 
gas phase. After evaporation, this molecule is believed to drive a fast high temperature 
gas-phase chemistry leading first to the formation of SO and SO2 and subsequently of other 
sulfur-bearing molecules such as CS and H 2 CS (Charnley 1997; Hatchell et al. 1998). It is 
still uncertain whether the release of H 2 S is caused exclusively by thermal evaporation or 
whether shocks play a role in the disruption of the grain mantles (Hatchell & Viti 2002). 
In low mass star-forming regions, where the temperatures are not high enough to evaporate 
the grain mantles, there is evidence that shocks may be responsible for injecting the H2S 
from ices into the gas-phase (Buckle & Fuller 2003). Sulfur-bearing molecules are considered 
to be an important tool for studying the presence of shocks within massive star-forming 
regions. Time-dependent chemical models have roughly succeeded to reproduce the observed 
abundances of the most common sulfur-bearing species both in high and low mass star- 
forming regions (Hatchell et al. 1998; Buckle & Fuller 2003). 

In starburst galaxies like NGC 253, only the most abundant sulfur-bearing molecules 
have been detected so far, namely CS, OCS, SO (Henkel & Bally 1985; Mauersberger et 
al. 1995; Petuchowski & Bennet 1992), and recently also S0 2 and NS (Martin et al. 2003). 
Extragalactic H 2 S has only been detected toward the Large Magellanic Cloud, where H 2 CS is 
also tentatively detected (Heikkila et al. 1999). In this paper we present the first detection of 
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H 2 S and the tentative detection of H 2 CS in NGC253, allowing us to investigate for the first 
time the sulfur chemistry in its nuclear region using a set of the most abundant sulfur-bearing 
species. 

2. Observations 

Most of the observations of the molecular species listed in Table 1 were carried out 
with the IRAM 30m telescope on Pico Veleta (Spain). The measurements were made in 
symmetrical wobbler switching mode with a frequency of 0.5 Hz and a beamthrow of 4' in 
azimuth. As spectrometers, we used two 256 x 4 MHz filterbanks for the transitions at 2 mm 
and two 512 x 1 MHz filterbanks for those at 3 mm. The J = 2 — 1 and J = 5 — 4 transitions of 
CS were measured with the 15 m Swedish-ESO Submillimetre Telescope (SEST) at La Silla 
(Chile). Observations were done in wide dual beam switch mode, with 11.8' beam throw 
in azimuth. Acousto Optical Spectrometers were used as backends with a resolution of 
1.4 MHz. Beam sizes ranged from 29" (at 85 GHz) to 12" (at 196 GHz) for the 30m telescope 
and from 51" (at 98 GHz) to 21" (at 244 GHz) for the SEST. In both telescopes, receivers 
were tuned to single side band with an image band rejection larger than lOdB. The spectra 
were calibrated with a dual load system. 

Fig. 1 shows the observed line profiles of CS and two of its isotopic substitutions ( 13 CS 
and C 34 S), H 2 S, SO, H 2 CS, and OCS, as well as the corresponding Gaussian fits. 

3. Results 

For most of the observed transitions, we identify two velocity components, which arise 
from molecular cloud complexes at two opposing sides of the nucleus (Mauersberger et al. 
1996), located at - 5" NE (the 180 km s" 1 component) and - 5" SW (the 280 km s" 1 com- 
ponent) as revealed by interferometric maps (Peng et al. 1996; Garcia-Burillo et al. 2000). 

The derived parameters from the Gaussian fits to the observed profiles are shown in 
Table 1. While most of the transitions of the main isotopic species of CS have already been 
reported (Mauersberger & Henkel 1989; Henkel et al. 1993), the improved signal-to-noise 
ratios of our data allow us for the first time to separate the two main velocity components 
in the line profiles. Unfortunately, the 2 mm lines of H 2 CS appear to be blended with other 
lines. The 4 14 — 3^3 transition is blended with the H36o; recombination line and the Ai 3 — 3i j2 
transition is contaminated by a contribution of the CS J = 3 — 2 line from the image band. 
In order to perform the Gaussian fits to the H 2 CS lines, we have respectively subtracted 
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a Gaussian profile similar to that observed for the H34a line measured at 160 GHz and 
the CS J = 3 — 2 emission measured in the signal band corrected by the corresponding 
image gain. As seen in the lower panel of the H 2 CS lines in Fig. 1, the CS and the H 2 CS 
lines do not fully account for the observed feature. The extra emission could be due to the 
presence of another unidentified line. In spite of the uncertainties due to the blending, both 
resulting H 2 CS profiles are self-consistent, since they are expected to have similar intensities. 
Nevertheless, this detection must be regarded as tentative. Similarly, the OCS J = 7 — 6 
transition appears blended with the J = 1 — line of HC 18 + . Additional observation of the 
HC 18 + J = 2 — 1 transition at 170 GHz allowed us to estimate and subtract the contribution 
of this species to the observed feature, by assuming local thermodynamic equilibrium (LTE) 
with an excitation temperature T cx ~12K. 



3.1. Molecular abundances and excitation temperatures 

In order to derive the physical properties of the molecular material we need to make 
an assumption on the total emission extent of the molecular emitting region. By smoothing 
the interferometric CS data from Peng et al. (1996), Mauersberger et al. (2003) derived the 
dependence of the observed line intensity on the beam size. This dependence indicates an 
equivalent source size of ~ 20" for the CS emitting region. The interferometer maps of CS 
only miss 36% of the flux (Peng et al. 1996) suggesting that the extent of the CS emission is 
not much larger than this estimate. As we have no additional information on the extent of 
the emission of the other rarer species, we will assume them to be similar to that of CS. This 
assumption is supported by the comparison of the high angular resolution maps of the CS 
emission with those of species with a fairly different chemistry such as SiO, H 13 CO + (Garcia- 
Burillo et al. 2000) and NH 3 (M. Lebron, private communication), where the morphology of 
the emission seems to be similar. 

We have derived source averaged column densities by correcting the measured line in- 
tensities for the estimated source size of 20" assuming optically thin emission. The corre- 
sponding population diagrams for CS, SO, and OCS are shown in Fig. 2. For CS we show 
the population diagrams for the two velocity components at 180 and 280 km s -1 and we 
derive a rotational temperature T rot ~ 10 K for both components. In these diagrams, the 
population in the J = 2 level, clearly above the linear fit, shows evidence for the presence 
of a lower excitation temperature component with T rot ~ 5 K which can be explained by 
the much larger region observed by the 51" beam of the SEST. As derived below, the CS 
emission is moderately optically thick. The effect of the opacity on the population diagram 
would result in a slight increase of the derived excitation temperatures. Thus we would get 
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temperatures of ~ 13 K and ~ 8 K for the two components. To derive the CS column density 
in Table 2 we have used the optically thin emission of one of its isotopes, 13 CS. For OCS we 
derive a T rot = 16 K. Our estimate of the OCS column density in Table 2 is in agreement 
with that derived by Mauersberger et al. (1995) when the differences in the assumed source 
size are taken into account. For the 180 km s™ 1 component of SO we obtain a T rot = 23 K. 
As seen in Fig. 2, the 280 km s -1 component of the SO 34 — 2 3 transition (E u /k ~ 29 K) 
has a higher intensity than expected if we assume, from the CS results, both components 
to have a similar rotational temperature (Fig. 2). This may be due to the contamimation 
of this component by the emission feature observed at a velocity of ~ 400 kms^ 1 which 
we tentatively identify as the 8i )8 — 7i i7 transition of NH 2 CN (Cyanamide). The measured 
integrated intensity has been then considered as an upper limit as shown in Fig. 2. We have 
therefore assumed a similar excitation temperature of 23 K for both velocity components to 
compute the column density of the 280kms _1 component of SO. Since we have only one 
transition for H2S and the two observed lines of H2CS have the same upper level energy, we 
need to make an assumption on the excitation temperature of these species. As the Einstein 
coefficients of these transitions are similar to those of CS and S0 2 , we will take an excitation 
temperature of 12 K for both species according to that derived from CS in this work, and 
from S0 2 by Martin et al. (2003). 

Table 2 shows the derived source averaged column densities, rotational temperatures and 
fractional abundances relative to H 2 for all the sulfur bearing molecules detected towards 
the nucleus of NGC253. The main sources of uncertainty in the derived abundances are 
those from the assumed extent of the emitting region and the H 2 column density, which are 
assumed to be the same for all the species. Assuming the extreme case in which the molecular 
emission is confined to a much smaller region of ~ 10" the derived source averaged column 
densities will be larger by less than a factor 2. However, these uncertainties are expected to be 
even smaller when column density ratios between different species are considered. A similar 
argument also applies to the uncertainties in the derived fractional abundances introduced 
by the assumed H 2 column density if, as already discussed, the spatial distribution among 
different species does not vary substantially. Since we will focus our discussion on the 
molecular abundance ratios between different species, these uncertainties will not affect our 
conclusions. 

Given that the J = 3 — 2 transition of CS, 13 CS, and C 34 S were observed with the 
same telescope we can compare the observed 12 CS to 13 CS line intensity ratio of 21 ± 3 
with the 12 C/ 13 C ratio of 40 ± 10 estimated by Henkel et al. (1993). Assuming that the 
12 CS/ 13 CS intensity ratio reflects the 12 C/ 13 C abundance ratio, we derive an optical depth 
for the 12 CS J = 3 — 2 transition of ~ 1.4. Therefore from the measured 13 CS/C 34 S line 
intensities and the 12 C/ 13 C ratio mentioned above we derive a 32 S/ 34 S abundance ratio of 
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8 ± 2. From our non detection of the C 33 S line, it is possible to give a 3a lower limit to the 
34 S/ 33 S ratio of 9. 

The line intensities corrected for beam dilution of the four observed transitions of CS 
and the J = 3 — 2 lines of 13 CS and C 34 S have been compared with the results of model 
calculations for the excitation of CS in which the Large Velocity Gradient approximation 
has been used. A linewidth of 100 km s -1 and a kinetic temperature Tki n =100K (similar to 
the rotational temperature of 100-142 K derived from NH 3 measurements by Mauersberger 
et al. (2003)) have been considered in these non-LTE calculations. Fig. 3 shows the results 
of this model for the observed line intensities and line ratios. The grey region corresponds 
to the best fit to the observed CS main isotope line intensities, calculated by minimizing the 
reduced \ 2 function as described by Nummelin et al. (2000). The best fit gives a CS column 
density N(CS)=5.7 x 10 13 cm~ 2 and a molecular hydrogen density n(H 2 )=2.2 x 10 5 cm" 3 for 
the velocity component at 180 km s" 1 and N(CS)=5.7 x 10 13 cm~ 2 and n(H 2 )=2.4 x 10 5 cm" 3 
for that at 280 km s -1 . These column densities are a factor of two lower than those in Table 2 
derived from the isotopomer 13 CS assuming LTE conditions. The H2 density derived from 
this analysis depends on the assumed kinetic temperature. Therefore, a T^ n of 50 K or 25 K 
would lead to H 2 densities a factor of 2 and 4 higher respectively. The choice of T^ n has, 
however, almost no effect on the column density determination. 

As seen in Fig. 3 shows, the best fit contour does not cover the J = 2 — 1 data at a 
probability level of 99.5%. This fact suggests that the CS emission is best explained if one 
consider several gas components in the region covered by the telescope beam. These compo- 
nents are related to the need of two different excitation temperatures to fit the population 
diagram. The J = 2 — 1 to 1 — ratio of 1.5 derived by Paglione et al. (1995), also plotted 
in Fig 3, supports this scenario. For the 180 km s -1 velocity component (Fig. 3(a)) a higher 
density component (n(H 2 )~ 4 x 10 5 cm -3 ) explains the J u = 5,4,3 transitions, and a lower 
density component (n(H 2 )~ 6 x 10 4 cm -3 ), with a 3 times higher column density, explains 
the lower J transitions. Similar results, within the errors, can be derived for the 280 km s -1 
component. 

From the non-LTE calculations we derive C 32 S/C 34 S ratios of 4-9 and 12 CS/ 13 CS ra- 
tios of 22-35 for the 180 and 280 kms" 1 velocity components respectively. The 12 CS/ 13 CS 
intensity ratio obtained for the 180 km s -1 component is half the 12 C/ 13 C ratio estimated by 
Henkel et al. (1993) while that obtained for the 280 km s -1 component is closer to their value. 
This may be explained by the fact that the 13 CS observation by Henkel et al. (1993) where 
made towards a position ~ 6" closer to the emission peak of the 280 kms -1 component. 
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4. Discussion 

4.1. Molecular fractional abundances 

We find that most of the sulfur-bearing species detected towards the nuclear region of 
NGC 253 show similar fractional abundances of a few 1CT 9 (Table 2). Only CS and OCS have 
abundances a factor of ~ 5 higher. In order to understand the mechanism that dominates 
the chemistry in the nuclear environment of NGC 253 we summarize in Table 3 the fractional 
abundances of all the sulfur-bearing molecules detected towards this starburst galaxy and 
those observed towards selected Galactic sources (hot cores, Galactic center clouds, dark 
clouds, photodissociation regions (PDRs) and shocks driven by massive protostars) which 
are dominated by different kinds of chemistry. The abundance ratios discussed on this section 
are also shown in Table 3. 

The SO2 abundance in PDRs is much lower,by a factor of 40, than that found in 
NGC 253. Thus, abundance ratios such as H2S/SO2 and CS/SO2 can be used to estab- 
lish the role of photodissociation in the chemistry of NGC 253. Particularly instructive is 
the SO/SO2 abundance ratio. SO/S0 2 ~ 100 towards PDRs is about two orders of mag- 
nitude larger than that in NGC 253 which suggests that photodissociation does not play a 
major role in the chemistry of its nuclear region. 

A comparison with dark clouds shows a clearly lower abundance of OCS than that 
observed towards NGC 253. The OCS/SO ratio is 0.1 and 0.4 in L134N and TMC-1 respec- 
tively, while we measure ~ 4 in NGC 253. The abundance ratios of S0 2 /SO and H 2 S/SO 
of ~ 0.2 and ~ 0.05 respectively measured towards these dark clouds are also much smaller 
than the corresponding ratios of ~ 1 observed towards NGC 253. Therefore, dark molecular 
clouds like L134N and TMC-1 do not seem to dominate the molecular composition near the 
nucleus of NGC 253. 

Chemistry driven by high velocity shocks in molecular outflows also differs significantly 
from the chemistry found in NGC 253. The OCS/SO ratio measured towards these sources 
is two orders of magnitude smaller than that in NGC 253, while CS/H 2 S is three orders of 
magnitude lower. The large abundance of H 2 S, as well as the low abundance of molecules 
such as CS, OCS, and H 2 CS in local molecular outflows, is not observed in NGC 253. 

As far as the hot cores are concerned, the molecular abundances vary significantly 
between Galactic sources. Different evolutionary states may account for these observed 
differences. As seen in Table 3 the abundance ratios in the Orion hot core mostly disagree 
with those in NGC 253. The relative abundances of most of the species towards SGRB2(N) 
are similar to what we find in NGC 253 except for the underabundance of H 2 S and OCS 
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with respect to the other molecules. The ratio OCS/SO is almost two orders of magnitude 
lower in both hot cores than in NGC253. In addition, the low rotational temperatures 
(T rot < 25K) derived for all the sulfur-bearing species in NGC 253 suggest that the observed 
emission cannot be due to hot cores where much higher rotational temperatures are expected 
(T kin > 70K). 

Martin et al. (2003) pointed out that the molecular abundances of SO2, NO, and NS 
are similar to those found in the envelope of the SGRB2 molecular complex. Table 3 also 
includes the comparison with two different clouds in the envelope of SGRB2. The first 
position, labeled SGRB2(OH), is located 30" south of SGRB2M, and the second one 100" 
north-east of SGRB2M. S0 2 , SO, and H 2 CS abundances in the Galactic center clouds agree 
well with those observed towards NGC 253. OCS/SO abundances ratios of 5 and 9 are similar 
to the ratio of ~ 4 found in NGC 253. 

Our study of the sulfur chemistry towards the nuclear region of NGC 253 seems to 
indicate a chemistry similar to that of the envelope of SGRB2 where low velocity shocks are 
thought to be the main heating mechanism of the molecular material (Flower et al. 1995; 
Mart in- Pint ado et al. 1997, 2001). The low rotational temperatures derived for NGC 253 
are also in agreement with those derived by Cummins et al. (1986) in the SGRB2 envelope. 
These results support the idea of large scale low velocity shocks driving the chemistry of 
the inner molecular material in the nucleus of NGC 253 which is consistent with the high 
rotational temperatures derived from NH 3 and the large observed SiO abundances (Garcia- 
Burillo et al. 2000; Mauersberger et al. 2003). 



4.2. Chemical models 

In order to get an insight into the dominant chemistry around the central region of 
NGC 253, we compare our results with those available from chemical models (Millar et al. 
1997; Charnley 1997; Hatchell et al. 1998; Buckle & Fuller 2003). For this we will assume that 
the whole nuclear region of NGC 253 can be described as a single molecular cloud. This is not 
likely as we are observing a 200 pc region where a number of unresolved molecular complexes 
with different physical and chemical properties coexist. Nevertheless, the dominant observed 
chemistry in the nuclear region of NGC 253 could be characterized by a single giant cloud 
as observed for the envelope of SGRB2 which has a size of ~ 30 pc. The SiO mapping of 
the Galactic center region show molecular cloud complexes extended over even larger scales 
(Martin-Pintado et al. 1997). 

We have compared the observed molecular abundances in NGC 253 with the time- 
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dependent models used by Hatchell et al. (1998) and Buckle & Fuller (2003) to describe 
the chemistry of hot cores and low mass protostars respectively. These models assume the 
release of H 2 S from grain mantles as the main precursor of the sulfur chemistry. Models with 
2kin < 50K and n(H 2 )< 10 5 cm -3 seem to best reproduce the abundances measured for SO2, 
SO, H2S, and H2CS shown in Table 3. This result is in agreement with the idea that the 
emission from the nuclear region of NGC253 is dominated by low excitation temperature 
and moderately dense molecular clouds similar to those found in the Galactic center clouds. 

However current chemical models fall short of accounting for the observed abundances 
of OCS by more than one order of magnitude. It has been pointed out by Hatchell et al. 
(1998) and van der Tak et al. (2003) that OCS may play an important role in the chemistry 
of sulfur-bearing molecules. Near-IR observations of massive protostars (Geballe et al. 1985; 
Tielens 1989; Palumbo et al. 1995) show that sulfur may freeze out onto grains not only in 
the form of H 2 S but also in the form of OCS. In fact OCS is the only sulfur bearing molecule 
detected so far in icy grain mantles. This would imply that OCS is also sputtered from 
grains by shocks. Hatchell et al. (1998) revised the chemical model by including an initial 
large abundance of OCS as well as reducing that of H2S, assuming these species are initially 
injected into the gas phase. These models, depending on the initially injected H 2 S and OCS 
(Wakelam et al. 2004), may account for the high OCS abundance observed in the nuclear 
region of NGC 253. 



4.3. Sulfur isotopic ratios 

The 32 S/ 34 S abundance ratios of 8 ± 2 and 13.5 ± 2.5 observed towards the nuclear 
starburst environment of NGC 253 (Sect. 3) and NGC 4945 (Wang et al. 2004) respectively 
are similar to those found within the Galactic inner 3 kpc region, but considerably lower than 
the value of ~ 24 measured in the Galactic disk (Chin et al. 1996). The estimated lower 
limit of the 34 S/ 33 S ratio towards NGC 253 is higher than the abundance of ~ 6 measured 
towards most of the Galactic sources studied by Chin et al. (1996). 

Massive stars, as well as type Ib/c and II supernovae, appear to slightly overproduce 
34 S and underproduce 33 S compared to 32 S (see Chin et al. 1996, and references therein). 
Thus, the compact stellar clusters observed near the center of NGC 253 (Watson et al. 1996) 
and the estimated overall star formation rate of S.QMqYt^ 1 (Strickland et al. 2004) and 
supernova rate (SNr) of 0.05-0.3 yr" 1 (Mattila & Meikle 2001; Ulvestad k Antonucci 1997) 
might account for the low 32 S/ 34 S as well as the enhanced 34 S/ 33 S ratios in NGC 253. This 
fact is supported by the high massive star formation rate of ~ 0.1M Q yr _1 estimated by 
Forbes et al. (1993) for the inner 6" of this starburst. 
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Uncertainties are large enough that the 32 S/ 34 S ratio observed towards NGC253 may 
not significantly differ from that in NGC4945, where slightly higher star formation and 
supernova rates (Strickland et al. 2004) might also have induced an overproduction of 34 S in 
its nuclear region. 
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Fig. 1.— Gaussian fits to the detected transitions of H 2 S, CS, 13 CS, C 34 S, SO, H 2 CS, and 
OCS towards the nuclear region of NGC 253 (a J20 oo = 00 h 47 m 33!4, <5 J20 oo = -25 17'23") The 
velocity resolution has been smoothed to lOkms -1 . Intensities are given in Tub (mK). The 
figure for H 2 CS shows in dashed lines the fits to the H 2 CS transitions, in dotted lines the 
assumed H36a and CS profiles and in thick line the total fit. Similarly, the figure for the 
OCS J = 7 — 6 transition shows the calculated HC 18 + in dotted line. 
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Fig. 2. — Population diagrams for the CS, SO, and OCS molecules. 180 km s -1 component 
is represented by squares and continuous lines, the 280 km s^ 1 component is shown with 
triangles and dashed lines. 
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Fig. 3. — LVG model predictions for the observed CS transitions and for the J = 3 — 2 
line of its isotopomers 13 CS and C 34 S at a T kin =100K as a function of the H 2 density and 
the CS column density. Upper and lower pannel correspond to the 180 and 280 kms" 1 
velocity components respectively. Lines show the measured line intensities and intensity 
ratios. The 2-1/1-0 ratio is from Paglione et al. (1995). Best fit to the line intensities of CS 
that minimizes the reduced x 2 function is shown in grey and corresponds to the probability 
of 99.5% of enclosing the true parameters. 
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a fixed parameter. 

b rms values in T m b and smoothed to a lOkrns -1 channel width. 
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Table 2: Physical parameters derived for each species. 
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Note. — Column densities are source averaged values over a 20"region. 
"Assumed N(H 2 ) = 1.710 22 cm~ 2 from Mauersberger et al. (2003). 

fc Column density derived using the optically thin 13 CS line and assuming a ratio 12 C/ 13 C=40 (Hcnkcl et al. 
1993). 

c from Martin et al. (2003). 
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Table 3: Relative abundances of S-bearing molecules in different environments. 
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Note. — Units of 10 9 are used. 
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(g) Minh et al. (1991); (h) Minh et al. (1990); (i) Martin et al. (2003); (j) McGonagle & Irvine (1997); (k) 
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